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A recently proposed statistical model for the effects of decoherence on electron transport manifests 
a decoherence-driven transition from quantum-coherent localized to ohmic behavior when applied 
to the one-dimensional Anderson model. Here we derive the resistivity in the ohmic case and 
show that the transition to localized behavior occurs when the coherence length surpasses a value 
which only depends on the second-order generalized Lyapunov exponent . We determine the 
exact value of of an infinite system for arbitrary uncorrelated disorder and electron energy. 
Likewise all higher even-order generalized Lyapunov exponents can be calculated, as exemplified for 
fourth order. An approximation for the localization length (inverse standard Lyapunov exponent) 
is presented, by assuming a log-normal limiting distribution for the dimensionless conductance T. 
This approximation works well in the limit of weak disorder, with the exception of the band edges 
and the band center. 

PACS numbers: 72.10.-d,72.15.Rn,71.23.An 



I. INTRODUCTION 

Anderson^ was the first to show that electronic eigen- 
states in disordered media can be localized, if the disor- 
der exceeds a certain threshold. In the infinite, discrete 
one-dimensional (ID) Anderson model 



(1) 



without correlations, where the onsite energies ei are in- 
dependently distributed according to a probability den- 
sity w{e) (with mean value 0, variance a^, third and 
fourth moments i^^ and k'^), all eigenstates are localized 
for cr^ > 0, i.e. the eigenfunctions decay exponentially in 
space. To simplify the notation we take t — 1 as energy 
unit and the lattice spacing a = 1 as length unit in the 
following. 

Even fifty years after its foundations, the physics of 
Anderson localization is a very active field of research, 
stimulated by the recent observations of localization of 
light and cold atoms, cf. Refs. 2-4 and references therein. 

The "absence of diffusion" due to disorder leads to an 
exponential increase of the electrical resistance of the 
localized system, if the electrons behave coherently.'^ 
Ohmic behavior, a resistance increase proportional to the 
linear extension of the device, can be achieved in spite of 
localization, if sufficient decoherence is included.EI In re- 
cent publications we have demonstrated this using a sta- 
tistical model for the effects of decoherence.'^'Sl Also the 
experimental observation of sequence-dependent conduc- 
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tance in DNA molecules (ohmic vs. exponential)"^ could 
be explained .121 

In this article we examine the decoherence-induced 
conductivity in the Anderson model subjected to the phe- 
nomenological model for decoherence. 

Decoherence-induced conductivity in disordered sys- 
tems was already subject of various papers. "'^^^^^ Here 
we examine how much decoherence is necessary in order 
to obtain finite conductivity and conclude that whenever 
the coherence length fulfills 



La < 



1 _ exp(-fi) 



(2) 



ohmic conductance is reached. Here is the second- 
order generalized Lyapunov exponent (GLE),'iSHi^ for 
which we report exact results for arbitrary diagonal dis- 
order and electron energy. All higher even-order GLEs, 
e.g. the fourth-order GLE x~^: can be calculated analyt- 
ically by the same method as The even-order GLEs 
have also been determined for nonlinear oscillators.'^^ 

Assuming a Gaussian limiting distribution for the log- 
arithm of the dimensionless conductance In T, we present 
an approximation in which the localization length A (the 
inverse standard Lyapunov exponent LE) is determined 
by ^ and x- This approximation works well in the limit 
of weak disorder with the exception of the band edges 
and the band center. 

The paper is organized as follows. Sect. [IT| explains the 
coherent transport formalism and the phenomenological 
statistical model by which we include decoherence effects. 
In Sect. |III| we derive the resistivity and the limitations 
under which finite conductivity is found. Sect. IV intro- 
duces the GLEs and derives Using the GLEs, in 
Sect. |V] we perform an approximate calculation of the 
localization length. Sect. |VI| discusses the results and 
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conclusions. The Appendices contain details of the cal- 
culations. 



II. MODEL 

The idea of the statistical model for the effects of de- 
coherence on electron transport is the following.'^ We de- 
scribe the electrons in a single-electron picture. Deco- 
herence is modeled by stochastic events: The electron, 
which in general is considered coherently, is extracted 
and re-injected at decoherence events, which take place 
at random locations in space. Thereby it loses its co- 
herence co mplete ly. The quantity of interest, e.g. the 
conductanc^'^HI] or resistanc^sH ig calculated for a given 
decoherence configuration^ a division of the system into 
coherent and decoherence regions. Afterwards the quan- 
tity of interest is averaged over an ensemble of decoher- 
ence configurations. Here we apply our statistical deco- 
herence approach to the Anderson model in the following 
way, cf. Fig. [T] which is a slight modification of the orig- 
inal. 

(a) 
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FIG. 1: (Color online) A coherent chain of A*' 
pled to wide-band contacts 



(a) 



sites cou- 
Three decoherence bonds 
introduced at random according to a probability p form a de- 
coherence configuration. |(b)| The resulting equivalent set of 
coherent chains. The denote the electron energy distribu- 
tion functions at the contacts and decoherence bonds. 



According to the probability p, the bonds of the An- 
derson Hamiltonian H are replaced by decoherence bonds, 
virtual contacts in the wide-band limit which cause a 
purely imaginary self-energy —irj to the neighboring re- 
gions of the system, forming a decoherence configuration. 
Thus p and rj are the two parameters by which we de- 
scribe phenomenologically the effect of decoherence. 

A decoherence bond represents the environmental de- 
cohering effect, or a decoherent escape into the rest of 
the electronic system. An electron is extracted out of the 
system at the decoherence bond (the bond represents a 
reservoir) and re-injected with the same energy (no in- 
elastic effects are considered) at the same bond involving 
a complete loss of phase coherence. 

The Green's function of a coherent subsystem which 
is delimited by two decoherence bonds and contains the 



sites reads 

G^[E~H + irj{\l){l\ + \j){j\)]-\ (3) 

where H is the Hamiltonian restricted to the sites between 
the two decoherence bonds, and E is the total electron 
energy. Due to the simple structure of the self-energies, 
the transmission through the subsystem according to the 
Nonequilibrium Green's function method (NEGF)^ reads 



V|(l|G|j)|^ 



(4) 



and can be evaluated recursively because of the tridiag- 
onal structure of G^^: 



1 _ |rj - V^Sj^i + i7?(rj_i + Sj) 



Tj 

with the polynomials 



(5) 



rj = {E- ej)rj^i - rj_2, Sj = {E - ej)sj-i - Sj_2, 
ro = 1, si = 1, 

r_i =0, So = 0. 

(6) 

For the derivation of Eqs. ([s]) and ([6]), see the Ap- 
pendix [Xj 

Although the dimensionless conductance T is often cal- 
culated without making reference to a coupling to exter- 
nal leads,!^ the results for large coherent systems which 
we discuss below do not depend on the model parame- 
ter 77, as it becomes an irrelevant boundary condition (cf. 
Appendix [C]) . 

Under the condition of a complete loss of phase coher- 
ence without energy relaxation at the decoherence bonds, 
we have shown in Ref. [71 that the resistance of a deco- 
herence configuration at infinitesimal bias is the sum of 
the individual subsystem resistances. 



h ^ 1 
2^4-1^' 



(7) 



g 2 

where k enumerates the coherent subsystems, and — = 
1 is the quantum of conductance assuming spin degener- 
acy, which we take as our conductance unit. To obtain 
the resistance as a system parameter, R is to be averaged 
over all decoherence configurations. 

Now consider the thermodynamic limit of an infinite 
system. As the locations of the decoherence bonds are 
uncorrelated, the relative frequency of the subsystem size 
j is pq^^^, where 



q = l-p 



(8) 



is the probability of a regular bond. On average, the 
subsystem size is 



I 

P 



(9) 



defining the coherence length in our model. In the infi- 
nite system, subsystems of size j appear infinitely many 
times, contributing on average {^/Tj) to the resistance. 
Here, (•) denotes the average over the disorder. 
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III. RESISTIVITY OF THE SYSTEM 

We define the resistivity of the system as the resistance 
per length in the hmit of an infinite system: 



i=i 



(10) 



Naturally this definition is restricted to those values 
q < q* for which the series converges. Then the system 
behaves ohmically. For q > q* , localization dominates, 
and the resistance increases exponentially with the sys- 
tem size. This is the decoherence induced conductivity: 
As the decoherence density p increases over the critical 
value p* — 1 ~ q* , there is a transition from localized to 
ohmic behavior. 

The disorder average {l/Tj) is a polynomial of order 

2j 



OO j 



-OO 

1 1 



(11) 



where we define 

°° 3 
Rj= Jde,...de,r]]lw{e,). (12) 

-00 »=i 

Using Eq. ^ we find a recursion for the Rj : 

00 

R, = Jde,... de, [{E' - 2Ee, + e^)r^_, 
-00 j 



= (£;2 + CT^) - 2ESj-i + Rj-2, 



where 



j 

^.^/d.,...de,.,.,_rn-(^0 



(13) 



(14) 



itself fulfills the recursion 



00 

Jde,...de,[{E^ 



ERj—i — Sj—i- 



(15) 



The initial conditions = 0, 5*0 = and Rq = I allow 
to determine the disorder-averaged resistance (l/Tj) of a 
subsystem of size j via Eqs. (13)-(15). 



Using the recursions in Eqs. ( 13 ) and ( 15 ) we calculate 



the generating function of the Rj 

00 

1 



(16) 



^1 



+ {E^-a^-l)z-E^- a^) 



where 

Ni = + il- E'^ + a^) z^ + (E'^ -1) z- 1. 



(17) 



With Eq. (Ill we evaluate the resistivity p [Eq. ( [10| ] 
in terms of the generating function TZ{z = q) arriving at 
the closed form valid for q < q* 



P 



2 

-p q 



2E^ 



4^2 4 



1 - (a^ + E'^) q-q^ 



47^2(1 + g) 
2^2 ^2 



(18) 



1 + g 



Fig. [2] displays the resistivity p as a function of the 
decoherence density p. 




FIG. 2: (Color online) The resistivity p as a function of the 
density of decoherence events p (=1/ L^, the inverse coherence 
length) for various disorder strengths a. For p < p* (p* are the 
poles in this graph) no finite conductivity exists. Parameters: 
Electron energy E — 0, broadening 77 = 1. 

The radius of convergence q* is given by the singularity 
of p which is nearest to the origin. I.e. q* is the smallest 
of the roots of A^i , which fulfill the relation 



1 2 



l-g2 



(1 + 9)^ 



tE' 



1. 



(19) 



We see that the singularities lie on an ellipse in the a-E 
plane, the semiaxes of which are determined by q. 
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Both summands of g{q) are strictly increasing for q G 
[0, 1), furthermore g{0) = and hniq_j.i g{q) = oo. Thus 
Eq. (19) obviously has a single, real solution q* £ (0, 1). 
The two other solutions are complex and fulfill |g| > 1, 
hence they have no physical meaning. This guarantees 
that for any disorder and any electron energy there is 
a critical decoherence density p* above which ohmic be- 
havior is achieved. 



The solutions qk of Eq. ( 19 1 read 
= -1) 



Qk 



3^ 



^ 1 

_ -2-^e-'i^-- [f2iE,a) + ^27h{E,a) 

where the functions 

fi{E,a) = -4: + 5E^ -E 
f2{E,a) = -16-24£;2 + 



(20) 



^2 , o 7712^2 4 

T + Z£j a — a 



\hE^ - 2E^ + 6(7^ - l2E^a^ 



ME,a) = 64E^ 



2E^a^ 



48E^ 



6E^a^ 
- 12E^ - 



4cr'* + 2QE^a'^ 



are polynomials and k — 0,1,2. Note that by ^- and 
(•)3 we denote the solution of = • and = • with the 
smallest complex argument arg(z) e [0, 2tt), respectively. 
Using ( 20 1 we can write 



q* = min{|qfc|} = 1 - p*. 
The dependence of p* on E is displayed in Fig. [3j 



(21) 
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FIG. 3: (Color online) The critical decoherence density p* for 
various values of cr as a function of the electron energy E. For 
p > p* , ohmic behavior is achieved, whereas for p < p* the 
resistance increases exponentially with the system size. In the 
limit (T — >■ we observe a kink at E = 2, the band edge of the 
system without disorder. 

For the interpretation of Fig. [3] note that the range 
Eg [—2,2] is the energy band of the system without dis- 
order. Thus only for \E\ < 2 there are eigenstates, which 



in this case are extended states. Electrons with \E\ > 2 
effectively have to tunnel through the system. This only 
leads to a finite resistivity if the tunneling length is re- 
duced by increasing the decoherence parameter p. The 
physical picture is similar for the case with disorder. Now 
there are no extended states, and electron transmission 
is exponentially suppressed for any E. Thus p* > for 
all values of E. 



IV. THE GENERALIZED LYAPUNOV 
EXPONENT 



The generalized Lyapunov exponents (GLEs) are de- 
fined aJ™ 



L(r) = - lim - In (t^ 

r i— >oo j \ ■' 



E.g. the second-order GLE £^ ^ = 2L{r ~ 2) reads 



= lim - In 



(22) 



(23) 



Since, according to Eq. (10 1, p/p^ is the generat- 
ing function of the (l/Tj), its pole q* determine^^Sl the 
asymptotic behavior of (l/Tj), and hence also 



lim 



r' = -Hq*)- (24) 



Therefore <^^^ can be determined exactly for any set of 
parameters {a,E) using Eqs. (20)~(21) and (24 1. 

Note that apart from the energy, ^"^ only depends on 
the second moment of the distribution of the onsite 
energies. Higher moments do not enter. 

Fig. [4] displays the dependence of on the disorder 
strength a for various energies E. 



0.0001 



I — I — 

5" =ln(£^/2-l+l0(£^/4 -1)'°) 




FIG. 4: (Color online) Dependence of ^ on a for different 
values of E. 

We observe that, depending on E and the disorder 
strength, assumes different functional dependencies 
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on a. In the weak disorder limit cr — ^ 0, ^ ^ approximates 

\E\<2, 

C'^{{2<y')\ \E\=2, (25) 

\E\ > 2, 



and for strong disorder cr — > oo 
as can be seen easily by inserting z 



(26) 

exp(— ^~^) into 
Eq. Jig). 

Outside the band \E\ > 2, is finite also for a = 0. 
This is why in Fig. [I] and Eq. ([25 1 no cr-dependence is 
given outside the band for weak disorder. 

Using we can restate the condition q < q* for a 
finite conductivity as: 



1 



(27) 



l-exp(-^i)' 

which for weak disorder inside the band becomes 

< C (28) 

As is a property of the infinite coherent system, 
it is determined purely by E and cr^ and does not de- 
pendon the decoherence parameters p and 77. Yet via 
Eq. (27), the second-order GLE determines the maxi- 



mum coherence length for which a finite conductivity can 
be reached. 



V. GENERALIZED LYAPUNOV EXPONENTS 
AND LOCALIZATION LENGTH 

The localization length A is defined as the typical 
length scale on which the eigenstates of a disordered sys- 
tem decay.'^B Its inverse (the standard Lyapunov expo- 
nent LE) is related to the dimensionless conductance of 
a coherent system via 



X-^ = - lim 



1 



(InT,). 



(29) 



It is easy to see that the localization length and the 
GLEs are connected by 



A" 



lim L(r). 

r— 



(30) 



With our method we can calculate L{r) only for even 
r 7^ 0, so that the limit r 0, Eq. ( [30| , is not obvious. 
It remains to be investigated in the future. 

We can, however, use our results to calculate A ap- 
proximately and gain new insight into the assumptions 
of single parameter scaling. 

One must distinguish two assumptions, which together 
define what is known as single parameter scaling (SPS) 
in localization theory.'^^^ISS 



Assumption A is that the variable Uj = — In Tj obeys 
the central limit theorem and approaches a Gaussian dis- 
tribution as j — > 00. 

Second, assumption B states that a single parameter 
suffices to characterize the Gaussian distribution, in other 
words, the mean Ci and the variance C2 are proportional 
to each other and that the proportionality factor is 2, 



C2/C1 = 2 for SPS. 



(31) 



Assumptions A and B were proven to hold for the con- 
tinuous Anderson model for uncorrelated weak disorder, 
provided complete pha se ran domization takes place be- 
tween scattering events. ^^^ISIl 

For the discrete Anderson model, however, already the 
weaker assumption A was shown to be wrong at the band 
edges and the band center.l^ Nevertheless, for weak dis- 
order the correct power law dependencies of A~^ on a 



are reproduced, if assumption A is applied, see Eqs. (44) 
and (45 ). Only the prefactors are wrong. For i?- values in 



between the band center and the band edges, assumption 
A even leads to the correct A^^ in leading order of tr, see 



Eq. (41) 



In the following we adopt assumption A as an approx- 
imation and investigate how strongly assumption B is 
violated. We find strong violation of assumption B at 
the band edges and at the band center. However, if for 
_E-values in between and weak disorder, one makes the 
assumption A, then assumption B is automatically ful- 
filled. 

The assumed Gaussian limiting distribution of the vari- 



able 



In Tj is characterized by its mean Ci and 



variance C2, where 



lim ^ = 2A-1. 



(32) 



This means that the probability P{uj) to find Uj € [u, u+ 
du] fulfills, for a finite system of length j, 

lim ( P{uj)duj e~' '^^2'' duj] ^ 0, (33) 

which allows to relate 
1 



^ = lim - In e" 

j^oo j 

1 



lim - In 



,Ci+C2/2 



(34) 



= lim - (Ci + C2/2) . 



where (•) denotes the average j_dujP{uj){-). 

Analogously to cf. Eq. ( |24[ ), we can also calculate 
the fourth-order GLE = 4L(4) (and similarly, all 
higher even-order GLEs), 




-ln(^*) 



(35) 



by determining z* , the pole of the generating function 
T-2{z) of (l/Tj) with the smallest absolute value. 



6 



The definition and derivation of 7'-2{z) can be found 
in the Appendix [B| Its denominator reads 



NixN2x{z^ 1), 



(36) 



where Ni is given in Eq. ( (TtI) an d N2 is a fifth-order 



polynomial in z given in Eq. (B29l 



Note that, as ^/T^ contains up to fourth powers of the 
polynomials Vj and Sj, the roots of N2, and therefore x 
itself, depend also on the third and fourth moments of 
the probability distribution w{e) of the onsite energies. 
At the same time, w{e) only enters through its moments 
into X, therefore any distributions with the same second, 
third and fourth moments has the same %, independently 
of the exact shape of the distribution. Analogously i(6) 
depends on the first six moments of w(e), etc. 

As A'2 is a fifth-order polynomial, the value of x can 
be calculated within numerical accuracy. Only for special 
cases (e.g. for E ^ 0, see below in Eq. (43 1), x can be 
represented in a closed form. 

Like we can relate x^^ to the cumulants C2 and 
Ci under the assumption A: 



X 



lim - In e^"j 
: lim - (2Ci + 2C2) . 



(37) 



This equation together with Eqs. (32) and (34) leads to 
the localization length 



under assumption A. 



(38) 



We solve Eqs. (34) and (37) for the ratio of the cumu- 
lants: 



C2 



2x-' - 4r 



4^ 



(39) 



We can compare this with the relation predicted by SPS 



theory, Eq. (31), by examining the parameter 



C2 



(40) 



which measures the relative deviation from SPS. Fig. [s] 
displays contour lines 6{(7,E) = const, for Gaussian on- 
site disorder, i.e. i/^ = 0, = 3a*. 

The interpretation of Fig. [5] is the following. For dis- 
order strengths a below the contour lines a < as (E) , the 
relative deviation from the SPS hypothesis is less than S. 
E.g. all points (a, E) below the black line in Fig. [5] fulfill 
5 < 10-4. 

For weak disorder and energies < \E \ < 2, we find 
X^^ = 3^-^ which according to Eq. (39) leads to SPS 
(Ca/Ci = 2). Using Eq. (l38|, we get 



8-2E^' 



for < |£;| < 2 and 0, (41) 



D 0.1 



0.01 




FIG. 5: (Color online) Contour lines for S = (C2/C1 - 2)/2 
which indicate the deviation from SPS. Assumed was a Gaus- 
sian onsite energy disorder. 



in agreement with the well-known formula by Thoules^^U 
which, in this energy range, is correct also for the discrete 
Anderson model.^ 

At the band center E — we find that SPS is violated 
even in the limit cr — > 0. Inserting 



= In 




for E = (42) 



as well as, for Gaussian onsite energy disorder k"* — Sct"*, 
3 



X 



In 



l + ^a^ 




for £; = 
(43) 



into Eq. (38) leads to the localization length 
2- v/3 



for i? = and a ^ 0. 



(44) 



This equation shows the correct power-law beh avior 
A^^ ^ a^ but deviates from the correct prefactor.'^^Iini 
This suggests that at the band center the assumption 
A breaks down and the simple approximation ( |38[ ) be- 
tween the localization length and the GLEs and x~^ 
is no longer valid. This is a signature of the band center 
anomaly.'23H2Sl 

Similarly, at the band edges \E\ — 2, SPS is violated 
also for cr — >■ 0. Eq. (38) leads to 



A" 



(2 



1/3 



1421/3) 



2/3 



for \E\ 



2 and a ^ 0, 
(45) 



again showing the correct power-law behavior A^^ ^ ct^/^ 
but deviating from the correct prefactor for the same 
reason as above. 
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Outside the band \E\ > 2, examining S does not lead 
to agreement with SPS for any disorder strength. 

In other words, if assumption A is apphed as an ap- 
proximation to the discrete Anderson model one obtains 
the correct localization length in leading order of a as 
long as E is in between the band edges and the band 
center. Then, assumption A automatically leads to ful- 
fillment of assumption B. 



VI. CONCLUSIONS 



rect weak disorder limit of the localization length. For 
this energy range, we present a maximum disorder as (E) 
which yields agreement with SPS, given a tolerance 6 of 
deviation. 

For the energies E — and = 2 (band edges and 
band center of the pure system), even for weak disorder 
the SPS hypothesis is not fulfilled. Furthermore, the lo- 
calization length one obtains under assum ption A does 
not agree with the literature results.l^^^ We conclude 
that at these energies assumption A breaks down, in ac- 
cordance with Ref. 1211 



The present article deals with the discrete-space ID 
Anderson model with uncorrelated diagonal disorder, 
characterized by its finite second through fourth mo- 



ments a 



and K*. 



When subjected to a statistical model for the ef- 
fects of decoherence (parameters p — the in- 
verse coherence length, and 77, the energy broadening 
due to decoherence) the length-dependent resistance 
at infinitesimal bias voltage undergoes a decoherence- 
induced transition from localized, exponentially increas- 
ing to ohmic, linearly increasing behavior, i.e. we find 
decoherence- induced conductivity. 

For the ohmic regime we derive the exact value of the 
resistivity p which for given model parameters p and 77 
only depends on the electron energy E and the second 
moment cr^ of the onsite energy disorder. Higher mo- 
ments do not enter. 

The critical decoherence density p* , at which the tran- 
sition from localized to ohmic behavior takes place, is 
connected to the second-order generalized Lyapunov ex- 
ponent (GLE) = -ln(l -p*). 

Therefore we can state the condition for finite con- 
ductivity as < (1 — exp(— ^~^))~^, which for weak 
disorder and energies inside the band becomes < ^, 
i.e. that the phase coherence length be smaller than the 
length ^. In this way, a property of the completely co- 
herent Anderson model (^) determines the limits [p*] of 
the ohmic behavior in the Anderson model including de- 
coherence. 

The GLE describes the exponential increase of 
(l/Tj), the disorder-averaged dimensionless resistance. 
Using the generating function of the {l/Tj) (which is 
nothing but p/p'^) we derive the exact value of for 
arbitrary electron energy and disorder. 

Analogously to the fourth-order GLE is calcu- 
lated from the generating function of (l/Tj^). In a similar 
way, all higher even-order GLEs can be determined. 

Under the assumption of a Gaussian limiting distribu- 
tion for InT (assumption A), ^ and x together allow to 
determine the localization length A and the ratio of the 
first two cumulants of InT, C2/C1 approximately. By 
studying this ratio we conclude that the single param- 
eter scahng (SPS) hypothesis which states C2/C1 = 2 
is automatically fulfilled in weak disorder for energies in 
between the band center and the band edges. Further- 
more in this energy range assumption A leads to the cor- 
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Appendix A: Derivation of the recursive formula for 
the Green function element Gij 

In this section we derive Eqs. (5]) and ^ which are 
necessary for the calculation of the transmission func- 
tion. The Green function of the finite chain with wide- 
band (virtual) contacts is the inverse of a tridiagonal, 
symmetric matrix: 



/E 



- £1 - 
-1 






i?7 -1 ••• 
£;-e2 -1 
-1 







-1 E-e, 



IT]/ 

(Al) 



As such, its components are easily derived by LU decom- 
position. 



/ai -1 

-1 '•• 


V : 

/I 

h 1 

'•• 



. -1 

1 ajj 

\ (mi -1 ...\ 

m2 



(A2) 







V; 1/^ 



V 



•• -1 

mjj 



(A3) 
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where a; = E-e^+iSnrj + iSiji] for i = 1, . . . , j, = - ^ 
for i = 1, . . . , j — 1, mi = ai and rrii = ai + li-i for 
*-2,...,j. 
From 



G-^G ^ LUG = 1 



(A4) 



follows that Gij is the first element of the vector x that 
fulfills 



L IJx^ — Sij . 

y 

Forward substitution yields 

y = Stj, 

and subsequent backward substitution in 

Ux = Sii 



gives 



Gij — xi — 



1 



nil X fn2 X • • • X nij 



(A5) 



(A6) 



(A7) 



(A8) 



rive a recursive expression for (1/Tj): 




16r;4 



oo 




J dei . . . dcj 




-oo 


.2 = 1 



1 

4 
1 



1 



2R 



2i?,_i + ii'R,-2) 
C.22 \ 



22 



where the i?, are defined in (12 1, and 



R',^lAr% 
Rf^iy^rU] 



— Ij + ll'^j+l^'j]; 
Imax{j,fc}[^jSfe] 

using the short-hand notation for the integration 



c22 



(Bl) 



(B2) 



(B3) 



Inserting the definitions of the m,j we find the recursive 
formula 



G 



1 



(A9) 



with the polynomials 

^0 = 1 

Pi=E -€i+ir^ 

P, = {E - (,)P,^i - P,^2, i-2,...,j-l 
P, = (i?-e,+ir;)P,_i-P,_2. 

As E, the Ei, and rj are real quantities, 



(AlO) 



(All) 



is the subdivision of the polynomial Pj into its real and 
imaginary parts, using the iterative polynomials and 
Si defined in (l6|. With the definition of the transmission 
function (Eq.^) follows Eq. 



Appendix B: Derivation of the generating function 

Using the relation ([s]) of the transmission function to 
the recursively defined polynomials rj and Sj (16]) we de- 



00 






jdei.. 


.dej 




—00 






Analogously we define 






R] = 







(B4) 



Ijb^'-i] 



(B5) 



and 



[iax{j,fc} ['"i'^i-l'^fc] 



112 



S'f " = lj[rjrj^iSjSj^i] 



(B6) 



By ^ the R and S fulfill the following recursions 

R!j — X^Rj_i — AX^Rj^i 



- AXiRf_i 



R 



■J -2 



6X2 i?f 

322 d13 

i-i 



rt — JLiJXj_i 



p22 



r?3i 



5: 



22 y q22 _\_ a V cllll I 022 

jj — ^4'Jj-ij-i -r '±^2»-'j-i 'r 'Jj-2j-2 



■211 



+ A:2(S'|iij_2 ^ 

- 2Xi(S'jly_2 

i22 



'S'j-2i-l) 
, c-211 



1) 



-1) 



9 



c22 




-1 


V c211 


-1 + Sj-lj- 


-2 


c211 

jj 


\^ o22 
- A35j_ij 


-1 




_ 1 + Xi Sj'^ 


2j-l 




Y 0211 


-1 




c211 

'^J-2j-l 




Cll2 


\^ o22 
- ^3^j-lj 


-1 




_1+Xi^22 


lj-2 




Y Cll2 






Cll2 

"^3-13-2 




Cr211 


\^ q22 


-1 


c211 






q112 


\^ q22 


-1 


Cll2 








\^ o22 
- A2ij_y 


-1 




, <^112 

-1 + '^j-lj- 





where we have defined 
E^ + 3£;f72 _ 



for n = 4 
for n = 3 
for n = 2 
for n — 1 



(B7) 

(B8) 
(B9) 



For the recursions, the following initial conditions apply: 



R\ 


— Xi, 


Rq - 




- -Ko 


n22 

- Rq - 


q22 


= ^2, 


Cll2 
"^11 


C22 


0211 

— 


q211 
— Dqi - 



C22 _ 1 
Oqi — 1 



(BIO) 



0112 _ rtllll 

•^10 ~ '-'1 



With this we can determine the generating functions of 
the R and S: 



= Xi + zXiU^iz) + Z + z^W^iz) 

- 4zX37^^l(z) + %zX-{R?'^{z) - 4zXl7^"(z) 

(Bll) 



7^3l(z) = ^i^fz^-l 

= X3 + zX37^''(z) - izX-{R?^{z) 
+ 3zXin^^{z) - zU^^iz) 



(B12) 



7^22(^)^^i^22^,-l 

= + zX27^''(z) - 2zXin^\z) + zU^^iz) 

(B13) 



(B14) 



= X2 + zXiSll{z) + 4zX25""(z) 

-2zX3(5fi"(z)+5ip(z)) 

+ zX2(522(z)+522(z)) 

- 2X,{Sl^{z) + Sll\z)) + z^Slliz) 

00 

= zX^Slliz) - 2zX^Sll^{z) + ^^^^(z) 



(B15) 



S 



S 



1111, 



(B16) 



Sll{z)^Y.S^ 



22 ^j-i 



1 + zX2Sll{z) - 2zXi5i2!i(z) + zSll{z) 



(B17) 



211/'^\ — 



zX35i\2(z) - 2zX2Sll''{z) + 



(B18) 



1127J-1 



: Xi + zX3522(z) - 2zX252ii(z) + zXiSlliz) 
zX2Sll\z) + 2zXi5""(z) - zSl^^iz) 

(B19) 



J = l 



.Z^-l 



= zXi5i\2(z)-z52ii(z) 

00 

5ii2(z)^E^-i-'"' 

= zXi5?2(z)-z5i\i2(z) 



(B20) 



(B21) 



(B22) 

: zX2SE{z) - + 5ii2(z)) 

z5iiii(z) 



Eqs. (B11|-(B14) and (B15)-(B22) constitute linear 



systems of equations for the generating functions TZ and 
S, respectively. 
The solutions read 

n\z) = [~z* + (-36x^X2 + 3X2 - Xi 

I\2 

+8X1X3 + 1 + - 6Xf + 2AXl) z^ 
+ (-6^2X4 + 3X2X4 - 8X1X3 + 12X^X2 
-3X2 - 18X2^ + 1 - AXl + X4 - 4X1 
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+24X1X2X3) + (-1 - 3X2X4 + 4X2 

+X4 - 6X2 _^ 4^2^ ^ _ ^323) 

7e22(z) = [(-2X2 ^ ^2 ^ ^2X1X3 - 3X2) ^ 

JV2 



and 



-Xo 



1 

' X2 
^Xa] 



(B24) 



[(-2X2 ^ ^2 ^ (2X1X3 - 3X2' 



(B25) 



1 



XiX2(z- 1) 



[-z^ + (24Xf - 36X2X2 



-4Xi2 + 8X1X3 + 6X| - X4 + 2X2 + 1) 
+ (-48Xi^ + 96Xf X2 + SXf - 16X1^X3 
-48Xi2x| - 4Xi2X4 + 4Xi2X2 + 32X1X3X3 
-12X| - 6Xi2 - 8X1X3 - 4X| + 2X2X4 
+2X| + X4 - 2X2 + 3) z^ + (32Xf + 12X^X4 
-72X^X2 - 48X1^X2X3 + 36X2 X^ _ 28^4 
+24X1^X3 + 8X2X1 - 12x2X2X4 + 24X2X1 
+24X1X22X3 - 18X| - 2Xi2X4 + 46Xi2X2 



12X1X2X3 



4X2X32 



-14Xi2 - 8X1X3 + 4Xf - 2X2X4 - 16X^ 
-2X4 - 4X2 - 3) + (-8X1^X4 + 8X^X2 
-32X5* X2X3 - 24X2X1 - 20Xf - 8X1^X3 



-12X,2X? 



10X12X2X4 + 16X^X22 



-4XiX|X3 + 3X^ + 6Xi2X4 - 2Xi2X2 
-28X1X2X3 + 4X2X1 - 3X|X4 + 14X| 
+8Xi2 + 8X1X3 + AXl - 2X2X4 - X22 



2X4 



+4X2 - 3) z'^ + (16X^ + 4X2X| - 2X2X2X 
-8X2x2 - 4X1X2X3 



3X, 



- 12X1 - 
3 



-HXfXa + 4X1X2X3 - 2X| 
-2X2X4 + 11X| - X4 + 2X2 + 3) z 
- (-2x2X4 + 2x2X2 + 4X1X2X3 - 2X1 

(2X2 



4X32 



-2X? - X? 



-X2 



l)z\ 

1 



X4 - 2X2 + 1) z2 



[-z^ + (20Xf - 32X2X2 



(B26) 

XiX2(z-l) L - ' v--^ —^-^ 
-8X1X3 + 5X22 - X4 + 2X2 + 1) z*^ + (-16Xf 



mxlxl - 4x2x4 

-3 



+48X1^X2 - 8X3X3 
+8Xi2X2 + 28X1X2X3 - 10X| - 4X1' 
-8X1X3 - 4X32 + 2X2X4 + 2X22 + X4 - 2X2 
+3) z^ + (4X1^X4 - 16X^^X2 - 16X1^X2X3 
+12Xi2x| - 16Xf + 8X1^X3 + 4X12X32 
-8X2X2X4 + l&XlXl + 20X1X^X3 - 15X^ 



+32x2X2 - 12X1X2X3 - 4X2X1 + 3X|X4 
+2X| + 12Xi2 - 8X1X3 + 4X| - 2X2X4 
-13X| + 2X4 - 4X2 - 3) z-* + {~AXf 
-4Xi2X| + 4X12X2X4 + 4Xi2x| - 4X1X22X3 
+3X^ + 4Xi2X4 - 4Xi2X2 - 20X1X2X3 
+4X2Xf - 3X22X4 + lOX^ + 4Xi2 + 8X1X3 
+4Xf - 2X2X4 - X22 - 2X4 + 4X2 - 3) z^ 
+ (-4X2X2 + 4X1X2X3 - 2X^ - 8X2 _ 4j^2 

+ 2X2X4 + 8X22 - X4 + 2X2 + 3) z2 + (-X2 

+X4 - 2X2 + 1) z - 1] (B27) 



where 



Xi 



(1 - 2X2 ^ ^2 _^ + X2) z - 1 



(B28) 



was already presented in Eq. (17) and 



7V2 = z^ + (-8X1X3 + 36Xi2X2 - 6X| - 24Xf 
-3X2 + 6Xi2 - 1 + X4) z'' + (6Xi2X4 + 3X2 
-3X2X4 - 2 - X4 + 18X| + 8X1X3 + 4Xi2 
-24X1X2X3 + 4X| - 12X12X2) z^ + (-10Xi2 



+6X32 - X4 + 2 - 4X32 + 3X2X4 + 3X2) 
+ (X4 - 3X2 + 1) z - 1. 



(B29) 



Here we have only displayed those generating functions 
which enter into the generating function of (l/Tj) defined 
below. The other functions TZ{z) and S{z) have similar 
expressions, particularly the same denominators. 



Using (B2 1 and the generating functions for the TZ and 
S (|B23)-(B27) we arrive at the generating function of 



5 1 1 



9 A 



+ 



4 1 - z V 4 16 
1 ^ -A ^4 

I6774 



+ 1^ + ^-) {n-\z)^sll{z)) 

(5i22(z)+5o2i2(z)), 



(B30) 



where 7?.(z) was defined in Eq. (16) 
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Appendix C: Calculation of the generating function Rjlj^i, IjRj^i . The solution reads 
of (l/Tj) without coupling to external leads 

Using the Green function of the tight-binding Ander- 
son Hamiltonian without couphng to the environment, 



T-i{z) = 



G= hm [E- H + irjlY 



(CI) 



1 being the unit matrix of adequate size, one defines the 
transmission probabihty of an electron between sites |1) 
and \j) as 



T,^\{l\G\j)\' 



(C2) 



Similar to Eqs. ^ and (|6| a recursive calculation of 
the transmission is possible: 



7V = 1^(^^+^)' 



where the Rj and Ij fulfill 



Rj — {E — ej)Rj_i — 77/j-i — Rj-2, 
Into the calculation of 



00 , 



(C3) 



(C4) 



(C5) 



enter the 7 generating functions of the disorder averages 
of the two-factor products R'j, RjRj^i, /|, Ijlj^i, Rjlj, 



-z^ + {E^ -a^- l)z -E^-a^ 

z3 + (cr2 -E^ + l)z2 + (^2 + E^ - l)z - 1 

(C6) 



and has the same poles as TZ in Eq. ( 16 ) 
Similarly, into the calculation of 



(C7) 



enter the 30 generating functions of the disorder averages 
of all four-factor products composed of Rj, Ij, Rj~i, and 

The closed form of T-2{z) is a very long expression, 
yet its denominator reads 



Nl X 7V| X (z - 1) 



(C8) 



with the same iVi and N2 as defined above in Eqs. (17) 
and (|B29|. 



Therefore, the generating functions T-i{z) and T-2{z) 
have the same poles q* and z* as reported above, and we 
find the same asymptotic behavior of (l/Tj) and (l/Tj) 
as j — 00. Hence the generalized Lyapunov exponents 
^"^ and of an infinite system are independent of 
whether one takes an environment into account or not. 
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